Transition magnetic moment of A — > S° in QCD sum rules 
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The A — > E° transition magnetic moment is computed in the QCD sum rules approach. Three 
independent tensor structures are derived in the external field method using generalized interpolating 
fields. They are analyzed together with the A and E° mass sum rules using a Monte-Carlo-based 
analysis, with attention to OPE convergence, ground-state dominance, and the role of the transitions 
in the intermediate states. Relations between sum rules for magnetic moments of A and E° and 
sum rules for transition magnetic moment of A — > E° are also examined. Our best prediction for the 
transition magnetic moment is ai e oa — 1.60±0.07 /ijv. A comparison is made with other calculations 
in the literature. 
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I. INTRODUCTION 

Magnetic moment is an important property in the elec- 
tromagnetic structure of baryons. In the case of octet 
baryons, various theoretical calculations have provided 
valuable insight into the underlying dynamics. The cal- 
culations benefited especially from the availability of pre- 
cise experimental information on the magnetic moments. 
The A — > E° transition magnetic moment is considered a 
member in this octet family, so any theoretical approach 
for octet baryon magnetic moments should be able to pro- 
vide an answer for this quantity. Relative to the attention 
paid to the individual members of the octet baryons, the 
A — > S° transition magnetic moment is less well-known. 

The QCD sum rule method is a nonperturbative ana- 
lytic formalism firmly entrenched in QCD with minimal 
modeling p|. In this method, hadron phenomenology is 
linked with the interactions of quarks and gluon through 
a few parameters: the QCD vacuum condensates, giving 
an unique perspective on how the properties of hadrons 
arise from nonperturbative interactions in the QCD vac- 
uum and how QCD works in this context. The method 
has been successfully applied to many aspects of strong- 
interaction physics. Calculations of the regular baryon 
magnetic moments have been carried out in QCD sum 
rules 0-Q ■ Limited study of the A — > S° transition mag- 
netic moment was performed in the traditional QCD sum 
rule approach [9( , and in the light-cone formulation of the 
approach 

In this work, we carry out a comprehensive, indepen- 
dent calculation of the transition magnetic moment of 
the A — > S° in the external field method in the tradi- 
tional QCD sum rule approach. This would complete 
the picture in the electromagnetic structure of the octet 
baryons from the perspective of QCD sum rules. A num- 
ber of improvements are made. First, we employ general- 
ized interpolating fields which allow us to use the optimal 
mixing of interpolating fields to achieve the best match 
in a sum rule. Second, we derive a new, complete set of 
QCD sum rules at all three tensor structures and analyze 
all of them. The previous sum rules, which were mostly 
limited to one of the tensor structures, correspond to a 



special case of the mixing in our sum rules. Third, we 
perform a Monte-Carlo analysis which can give realis- 
tic error bars and which has become standard nowadays. 
Fourth, we analyze the three sum rules using the A and 
E° mass sum rules derived from the same generalized in- 
terpolating fields. The performance of each of the sum 
rules is examined using the criteria of OPE convergence 
and ground-state dominance, along with the role of the 
transitions in intermediate states. Finally, we examine 
the relations between magnetic moment of A and S° and 
transition magnetic moment of A — > E°, using the newly- 
derived QCD sum rules. 



II. METHOD 

The starting point is the time-ordered correlation func- 
tion in the QCD vacuum in the presence of a static back- 
ground electromagnetic field F^: 



n(p) = i / d 4 x e lp - x (0 | T{ if (x) r? A (0) } | 0)j 



(1) 



On the quark level, it describes a hadron as quarks and 
gluons interacting in the QCD vacuum. On the phe- 
nomenological level, it is saturated by a tower of hadronic 
intermediate states with the same quantum numbers. By 
matching the two descriptions, a connection can be es- 
tablished between hadronic properties and the under- 
lying quark and gluon degrees of freedom governed by 
QCD. Here r\ is the interpolating field that couples to 
the hadron under consideration. The subscript F means 
that the correlation function is to be evaluated with an 
electromagnetic interaction term added to the QCD La- 
grangian: Ci 
tromagnetic potential and 
tromagnetic current. 

Since the external field can be made arbitrarily small, 
one can expand the correlation function 



-A^J^, where is the external elec- 
eq(p{^q is the quark elec- 



n(p) = ii^{p) + n < - 1 \ P ) 



(2) 



where Iv- \p) is the correlation function in the absence 
of the field, and gives rise to the mass sum rules of the 
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baryons. The transition magnetic moment will be ex- 
tracted from the QCD sum rules obtained from the linear 
response function U^(p). 

The action of the external electromagnetic field is two- 
fold: it couples directly to the quarks in the baryon in- 
terpolating fields, and it also polarizes the QCD vacuum. 
The latter can be described by introducing new parame- 
ters called vacuum susceptibilities. 

The interpolating field is constructed from quark fields 
with the quantum number of baryon under consideration 
and it is not unique. We consider a linear combination 
of the two standard local interpolating fields. 

v A (uds) = ~2^le abc [2(u aT Cj 5 d b )s c + {u aT C l5 s b )d c 

-{d aT Cj 5 s b )u c + P(2(u aT Cd b ) lb s c 
+(u aT Cs b )~f 5 d c - (d aT Cs b )j 5 u c )}; 

(3) 



r] sa (uds) = -V2e abc [{u aT Cj 5 s b )d c + (d aT C-/ 5 s b )u c 
+P{(u aT Cs b ) l5 d c + {d aT Cs b ) l5 u c )]. 

(4) 

Here C is the charge conjugation operator; the super- 
script T means transpose; and t a bc makes it color-singlet. 
The normalization factors are chosen so that correlation 
functions of these interpolating fields coincide with each 
other under SU(3)-flavor symmetry. The real parameter 
(3 allows for the mixture of the two independent currents. 
The choice advocated by Ioffe [2| and often used in QCD 
sum rules studies corresponds to (3 — —1. We will take 
advantage of this freedom to achieve optimal matching 
in the sum rule analysis. 



I 

A. Phenomenological Representation 



We start with the structure of the two-point correlation function in the presence of the electromagnetic vertex to 
first order 

U(p) = i J d 4 xe^ < 0|ry E °(x)H / d 4 yA^y)J»(y)}f} A (0)\0 > F . (5) 

Inserting two complete sets of physical intermediate states, we restrict our attention only to the positive energy ones 
and write 

n(p) = / d A xd A yj^ (|^4 E W -m\, -ie -Ml-ie , g s 

e ipx A^(y) < Ol^'ix^N'k's' >< N'k' s'\J»(y)\Nks >< Nks\fj A (0)\0 > . 
The matrix element of the electromagnetic current has the general form 

<fcV|J^O)|fc.s>=^ S 0[^i^^^ (7) 

where q — k' — k is the momentum transfer and Q 2 = —q 2 . In the following, we treat A and E° mass as degenerate 
because their mass difference is small and set m = mA + m s° to be the average mass of A and S°. The Dirac form 
factors Fx and F 2 are related to the Sachs form factors by 

G £ (Q 2 )=F 1 (Q 2 ) + ^F 2 (g 2 ) 
G M (Q 2 )^F 1 (Q 2 )+F 2 (Q 2 ). 

At Q 2 = 0, -Fi(O) = 1, ^2(0) = /i^oA w hich is the anomalous transition magnetic moment, and Gm(0) = -Pi(O) + 
-^2(0) = Ms°A which is the transition magnetic moment. Note that we consider here the transition A — > £ . The 
moment for the reverse transition is related by a minus sign: /iAS° = — 

Writing out explicitly only the contribution of the ground-state nucleon and denoting the excited state contribution 
by ESC, we arrive at 



n(p) - |A A A s oF^[p 2 - m 2 - ie]- 2 (p + m){ ^|^ a^(p + m) + io*» 
-(p M 7" -P v l ll )(p + m)[p 2 - fh 2 - is}- 1 } + ESC, 



(9) 



where Aa and A^o are phenomenological parameters that measure the overlap of the interpolating fields with the ground 
state. Examination of its tensor structure using Mathematica reveals that its has 3 independent combinations, which 
we name with the following short-hand notation: 



WEx = F^ipa^ + a^p), WOi = F^a^, W0 2 = F""i(p^ - PuT»)P- 



(10) 
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They are the same as those for the magnetic moment of octet baryons Q . Upon Borel transform the ground state 
takes the form 



4M Z m 



2(/"E»A - 1) 



[W0 2 ]}, (11) 



where M is the Borel mass. 

Here we must treat the excited states with care. For a generic invariant function, the pole structure in momentum 
space can be written as 



°iv<->jy _(_ n ^ N* 

(p 2 ~fh 2 ) 2 j-^ (p 2 -m 2 )(p 2 -m* 2 ) 



E 



N 



J (p 2 —rn* ) 



(12) 



where Cn^n, CWoiV* and Cn*^n* are constants. The first term is the ground state pole which contains the desired 
transition magnetic moment The second term represents the non-diagonal transitions between the ground state 

and the excited states caused by the external field. The third term is pure excited state contributions. Upon Borel 
transform, it takes the form 



AaA s qm s q a -m 2 /M 2 1 „-m 2 /M 2 
M 2 e t c 



E 



n 2 \ 



1 



E 



M 2 



2 /M 2 



(13) 



The important point is that the transitions give rise to a contribution that is not exponentially suppressed relative to 
the ground state. This is a general feature of the external-field technique. The strength of such transitions at each 
structure is a priori unknown and is an additional source of contamination in the determination of the transition 
magnetic moment /is A- The standard treatment of the transitions is to approximate the quantity in the square 
brackets by a constant, which is to be extracted from the sum rule along with the ground state property of interest. 
Inclusion of such contributions is necessary for the correct extraction of the magnetic moments. The pure excited 
state contributions arc exponentially suppressed relative to the ground state and can be modeled in the usual way by 
introducing a continuum model and threshold parameter. 



B. Calculation of the QCD Side 

On the quark level, by contracting out the quark pairs in Eq. ([1]) using the interpolating fields in Eq. ([3]) and 
Eq. (U), we obtain the following master formula in terms of quark propagators, 

(n I T{ if {x) 7? A (0) } I ft) = -2e abc e a ' b ' c '{ 

~sf l5 CSf T C l5 sf + 2Sf l5 CS h f C^Sf' - 2Sf l5 CSf T C l5 sf 

+Sf l5 CSf T C l5 sf + sfTr[ l5 CSf T C l5 S c u c '} - Sf Tr[ lb CS bb ' T C l5 S c d c '\ 

-pj 5 sf l5 CSf T CSf ~ (3S a /CS b s b ' T C l5 S c d c ' l5 + 2(3 l5 sf l5 CS b u b ' T CSf 

+2psfCS b u b ' T C l5 Sf l5 - 2(3 l5 sf l5 CS bb ' T CSf - 2(3Sf 'CS b f C l5 Sf 75 

-(3 2 l5 sfCS bb ' T CS c / l5 + 2f3 2 l5 sfCS b u b ' T CS c / l5 - 2p 2 l5 sf CS b f CS c u c ' l5 

+(3j 5 Sf l5 CS bb ' T CSf + pSfCS b f 'c l5 Sf l5 + (3 l5 SfTrl l5 CS bb ' T CSf] 

+(3Sf l5 Tr[CS bb ' T C l5 Sf] - ' Tr[ l5 C S b / ' CSf] - flS? ' l5 Tr[C S b / ' C l5 S c /} 

+/3 2 l5 S a /CS bb ' T CSf l5 + p 2 l5 Sf l5 Tr[CS bb ' T CSf) - f5 2 l5 S a / l5 Tr[CS bb ' T CSf], (14) 

where S^ b (x, 0; F) = (0 | T{ q a (x) q b (0) } | 0)f is the fully interacting quark propagator in the presence of the electro- 
magnetic field. The propagator can be derived from the operator product expansion (OPE). To first order in F^ v and 
m q (assume m u — nid = 0, m s 7^ 0), and up to x 4 , its expression can be found in Ref. |8j. 



C. The QCD Sum Rules 



With the above elements in hand, it is straightforward to evaluate the correlation function by substituting the quark 
propagator into the master formula. We keep terms to first order in the external field and in the strange quark mass. 
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Terms up to dimension 8 are considered. The algebra is very tedious. Each term in the master formula is a product 
of three copies of the quark propagator. There are hundreds of such terms summed over various color permutations. 
We used a Mathematica package MathQCDSR [ll[ that we developed to carry out the calculations. We confirmed 
that the QCD side has the same tensor structure as the phenomenological side. 

Once we have the QCD side (LHS) and the phenomenological side (RHS), we can construct the sum rules by 
matching the two sides in the standard way. Since there are three independent tensor structures, three sum rules can 
be constructed. 

At the structure WEi, the sum rule can be expressed in the following form 



Cl L- 4 / 9 E 2 (w)M 4 + c 2 m sX aL~ 2 ^ 27 E 1 (w)M 2 + c 3X a 2 L^ 27 E (w) + c 4 bL- 4 / 9 E (w) + (c s + c 6 )m s aL- 4 / 9 E (w) 
+(c 7 + c s )a 2 L 4 / 9 ^ + Cg xm 2 a 2 L-^ 2 ^ + c w m s m 2 aL- 2 ^ 27 'ik 

(15) 

where the coefficients are given by: 



Cl = i(2(/3 2 + /3 + l)e u - ((3(/3 + 2) + 2)e d ) 

c 2 = |/?(3(/3 + l)e u - (2/3 + 3)e d ) 

c 3 = ~(l)((3-l)(e d -e u )((3 + f s + l) 

c 4 = i((e s + 5e u )f3 2 + 4e u f3 - (J3 + 2) 2 e d + 4e„) 

c 5 = Me d (~6f s + Pifi + 308 + 4)/ s + 10) - 2) - 2e„(-3/ s + /3 (f3 + l)(6/ s + 5) - 1)) 

c 6 = MSeJMP 2 + ((3 2 +f3- 2)e d (2 K - + e u ((((3 - 2)f3 + 4)k + ((3 2 + f3 - 2)0) (16) 

ct = Ts(-e d ((3 2 - 2f s (3 + 5/ s - 1) - (fi - l)e u (5/ s + 0(3 /. - 1) - 1)) 

c 8 = i(3e s / s /^/3 2 + {fi- l)e d ((5f s + (3(3 f. - 1) - 1) « - (fi + f. + 1)0 

+e u (((3 2 - 2f s (3 + 5f s - 1)k + (fi - 1) ((3 + f s + 1)0) 
c 9 = ^(esfsW 2 + ((3- l)e d (5f s + 2{3(f s + 2) + 4) - e u ((f s + 4) f3 2 + 3/ s /3 - 5/ s - 4)) 
cio = J((5e»(/ S + 1) - 2e s )/3 2 + 6e„(/ s + l)/3 - (/3(/3 + 6) - 12) e d (f s + 1) - 12e u (/ s + 1)). 



In addition to the standard quark condensate, gluon condensate, and the mixed condensate which are rescaled as 



a = -(2tt) 2 (mi), b = (g 2 G 2 ), (ug c a ■ Gu) 



Three vacuum susceptibilities caused by the external field 



are introduced: (qa f _ lu q) F = e q x{qq)F llv , (qg c G t _ LU q) F = e q K{qq)F llu , (qg c e^ P xG pX "f 5 q) F = ie q ^{qq)F fll/ . Note that 
X has the dimension of GeV~ 2 , while k and £ are dimensionless. The quark charge factors e q are given in units 
of electric charge: e u = 2/3, e d = —1/3, e s = —1/3. Note that we choose to keep the quark charge factors 
explicit. The advantage is that it can facilitate the study of individual quark contribution to the transition magnetic 
moment. The parameters / and <fr account for the flavor-symmetry breaking of the strange quark in the condensates 
and susceptibilities: f = -kA = ^ 9 <= cr '^ s \ A, = K± — _ ^ ne arLO malous dimension corrections of the 
interpolating fields and the various operators are taken into account in the leading logarithmic approximation via the 



factor L 1 



a s (M 2 ) 



\n(M 2 /A 2 QCD ) 

M p 2 /k 2 qcd ) 



where fi = 500 McV is the renormalization scale and Aqcd is the QCD 

scale parameter. As usual, the pure excited state contributions are modeled using terms on the OPE side surviving 

M 2 — > oo under the assumption of duality, and are represented by the factors E n (w) = 1 — e~ w l M ^ > 
where w is an effective continuum threshold and it is in principle different for different sum rules and we will treat 
it as a free parameter in the analysis. Also, A is the rescaled current coupling A = (27r) 2 A. The rescaling is done so 
there is no explicit factors of ir in the sum rules. 



At the structure WOi, the sum rule can be expressed in the following form 



Cl m s L- 8 / 9 E 2 (w)M 4 + c 2X aL- 16 / 27 E 2 (w)M 4 + (c3+c 4 )a£i(u>)M 2 + 
+c 6 xabL- 16 / 27 E (w) + (c 7 + Cg)™^ 2 ^ + c 9 ab^ 
= -A A A s om r^JSA 4- ^aHAZi 4- A~\ p-^M 



c 5 m s xa 2 L- 1 ^ 27 E (w) 



M 2 



(17) 
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where the coefficients are 

ci = (P 2 + p-2)(e d -e u ) 

C2 = i((4/3 2 +p- 5)e u - {/3 2 + p - 5)e d ) 

c 3 = |(e d (-3/ s + /3(/3(3/ s - 2) + 4) + 1) - (p - l)e M (-5/3 + 3 (/? + l)/ s - 1)) 

c 4 - i(G8 - l)e B (2(8 ) 8K + k) + (0 + 8)0 + e d (2((7 - 5/3) + 1)« + (£(2/3 - 7) + 8)0) 

c 5 = |((/3 2 + /3 + l)e d - (2[3 2 + p + l)e„)(/ s + 1) (18) 

c 6 - Tsll" 2 / 32 + P + 4 ) e <* + (Z 3 " + 4 K) 

c 7 = A(e d ((7/ S + 4)/3 2 + 4/ s /3 + /? - 2/ s + 4) - e„((10/ s + 7) /3 2 + 4/ s /3 + /3 - 2/ s + 4)) 
c 8 = i(e d (/ s + P(P + 4(0 + l)/ s + 4) + 4) - e u (/ s + /3(7/ s /3 + 4 /3 + 4/ s + 4) + 4))(2k - 
c s = 558 ((- 2 ^ 2 + P + !) e d " e " " ~ e «) + e «))- 

At the structure WO2, the sum rule can be expressed in the following form 

Cim s L- 8 / 9 £;i(w)M 2 + c 2 xaL- 16 / 27 E 1 (w)M 2 + {c 3 + Ci )aE a {w) + c 5 m sX a 2 L- 16 ^ 7 E (w) + c^aL-^ 9 
+c 7 xabL- 16 ^ 7 + (c 8 + c 9 )m s a 2 ^+c w ab^ (19) 

_ AaA s q r 2(^ s o A -l) A]„-m 2 /M 2 

- rfT~ i rAje 

where the coefficients for A — > S° at WO2: 

ci = i(/? 2 +/?-2)(e d -e„) 
c 2 = f(/3-l) 2 (e d -e u ) 

cs = |M/? 2 (3/ s + 2) - 3(/ s + 1)) - 3(/? 2 - l)e u (f s + 1)) 

c 4 = i(3(/3 - l)e„(/?£ - 2k) + e d ((3 - 2/3)/3£ - 2 ((/? - 3)/? + 3)/c)) 

c 5 - |(/3 - l)(e rf - e u )(-/3 + (/3 + 3)/ s + 3) 

c 6 = s((-ea + e u + 4e„/ s )/3 2 + 3e u f3 - e u (4f s + 5)-(p- l)e d (2 p + 4(j9 + l)/ s + 5)) lZUj 

= T | 3 (/3-l) 2 (e (i -e u ) 
c 8 - i(e u ((10/ s + 7)/3 2 + Af s p + p-2f s + 4)- e d ((7f s + 4) p 2 + 4f s p + p-2f s + 4)) 
c 9 = -(i)(e d (/ s + p(P + 4{p + l)f s + 4) + 4) - e u (fs + P(7f s P + 4 P + 4/ s + 4) + 4))(2/s - 
cio = 25s ((/3 - l)(2/3 + l)e d + e u + P(P(e s - e u ) + e„)). 

Note that the sum rule from WEi involves only dimension-even condensates, so we call this sum rule chiral-even. 
The sum rule from both WOi and WO2 involves only dime nsion-odd condensates, so we call them chiral-odd. 



III. SUM RULE ANALYSIS 

The sum rules for transition magnetic moment have 
the generic form of OPE - ESC — Pole + Transition, or 

n A ^ E c (QCD, P, w, M 2 ) = A A A s o + A) e-^' M * , 

(21) 

where QCD represents all the QCD input parameters. 
The mathematical task then boils down to the following: 
given the function n^x: with known QCD input pa- 
rameters and the ability to vary /?, find the phenomeno- 
logical parameters (mass m, transition magnetic moment 
/is A) transition strength A, coupling strength Aa and 
As" , and continuum threshold w) by matching the two 
sides over some region in the Borel mass M . A \ 2 mini- 
mization is best suited for this purpose. It turns out that 
there are too many fit parameters for this procedure to 
be successful in general. To alleviate the situation, we 
employ the corresponding mass sum rules for the parti- 
cles [1, [lH which have a similar generic form of OPE - 
ESC = Pole, or 

n A (QCD,P,w u M 2 ) = \\e- m ^ M \ (22) 



IIj:o(QCD,P,w 2 ,M 2 ) = \\ a e- m ^l M \ (23) 

They share some of the common parameters and factors 
with the transition sum rules. Note that the continuum 
thresholds may not be the same in different sum rules. 
By taking the following combination of the transition and 
mass sum rules, 

U A ^o(QCD,P,w,M 2 ) = fi^A | 

y/U A (QCD, p, wi , A/ 2 )n s o (QCD, p, w 2 ,M 2 ) M 2 

(24) 

the couplings A and the exponential factors are canceled 
out. This is the form we are going to implement. By plot- 
ting the two sides as a function of 1/M 2 , the slope will 
be the transition magnetic moment ^s°A an d the inter- 
cept the transition strengthA The linearity (or deviation 
from it) of the left-hand side gives an indication of OPE 
convergence and the role of excited states. The two sides 
are expected to match for a good sum rule. This way 
of matching the sum rules has two advantages. First, 
the slope, which is the transition magnetic moment of 
interest, is usually better determined than the intercept. 
Second, by allowing the possibility of different continuum 
thresholds, we ensure that both sum rules stay in their 
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valid regimes. 

We use the Monte-Carlo procedure @, [H, [l4| to carry 
out the search. In this method, the entire phase-space 
of the input QCD parameters is explored simultaneously, 
and is mapped into uncertainties in the phenomenologi- 
cal parameters. This leads to more realistic uncertainty 
estimates than traditional approaches. 

The QCD input parameters are given as follows. The 
condensates are taken as a = 0.52 GeV 3 , b = 1.2 GeV 4 , 
rn,Q = 0.72 GeV 2 . For the factorization violation pa- 
rameter, we use k v — 2.0. The QCD scale parameter is 
restricted to Aq C d = 0.15 GeV. The vacuum suscepti- 
bilities have been estimated in studies of nucleon mag- 
netic moments @, H, [TH, [l6[ , but the values vary depend- 
ing on the method used. We use x = —6.0 GeV~ 2 and 
K = 0.75, £ = —1.5. Note that x is almost an order of 
magnitude larger than k and £, and is the most important 
of the three. The strange quark parameters are placed 
at m s = 0.15 GeV, f = 0.83, = 0.60 0, H|. These 
input parameters are just central values. We will explore 
sensitivity to these parameters by assigning uncertainties 
to them in the Monte-Carlo analysis. 



IV. RESULTS AND DISCUSSION 
A. Transition magnetic moments of A — ¥ S° 

We analyzed all three sum rules. For each sum rule, 
we have in principle 6 parameters to determine: //£°Aj 
A, w, u>i, u>2, /3. But a search treating all six parameters 
as free does not work because there is not enough infor- 
mation in the OPE to resolve them. Here, we use the 
freedom to vary /3 as an advantage to find the optimal 
match. Another parameter that can be used to our ad- 
vantage is the continuum threshold w\ for the mass sum 
rule. We fix it to the value that gives the best solution 
to the mass sum rule independently. For A, w\ = 1.60 
GeV; for E°, W2 — 1-66 GeV. This way the transition 
magnetic moment sum rule and the mass sum rule can 
stay in their respective valid Borel regimes. This leaves 
us with three parameters: /Lts°Ai A, w. Unfortunately, a 
three-parameter search is either unstable or returns val- 
ues for w smaller than the particle mass, an unphysical 
situation. Again we think this is a symptom of insuf- 
ficient information in the OPE. So we are forced to fix 
the continuum threshold w that corresponds to the best 
match for the central values of the QCD parameters. 

The Borel window is determined by the following two 
criteria: OPE convergence and ground-state dominance. 
It is done iteratively. For each value of (3, we adjust the 
Borel window until the best solution is found. Based on 
the quality of the match, the broadness of the Borel win- 
dow and its reach into the lower end (non-perturbativc 
regime), the size of the continuum contribution, and the 
OPE convergence, we find that the chiral-odd sum rule at 
structure WO2 is more reliable than the chiral-even sum 
rule at structure WEi and chiral-odd sum rule at struc- 



ture WOi. The physical reason is that the contributions 
of positive- and negative-parity excited states partially 
cancel each other in the chiral-odd sum rules whereas 
they add up in the chiral-even sum rules. Since the WO2 
sum rule has power corrections up to 1/Af 4 while WOi 
sum rule power corrections only goes up to 1/Af 2 , the 
WO2 sum rule is expected to be more reliable than the 
WOi sum rule. In the following, we present some de- 
tailed analysis of this sum rule in order to demonstrate 
how it produces the results. 

Fig. [U shows the matching of the two sides in Eq. (|24l) . 
The left hand side is normalized OPE minus ESC con- 
tribution and the right hand side is pole plus transition 
contribution. According to the right-hand side of this 
equation, the plot should be linear as a function of 1/Af 2 , 
so we try to match the left hand side over a Borel window 
by a linear line. We can see that the two curves agree 
very well over the Borel window 1.2 to 1.7. Based on 
this match, the slope 1.60 gives the transition magnetic 
moment /i£0A, and the intercept 0.17 gives the transition 
contribution A between intermediate states. This match 
was performed at j3 = —0.2 and w = 1.70 GeV. 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

M~ 2 (GeV -2 ) 



FIG. 1: The matching for the chiral-odd QCD sum rule at 
structure WO2 in Eq. (TE?)) as prescribed in Eq. (J2U). The 
error band is derived from about 2000 QCD parameter sets 
in the Monte Carlo analysis. 

Figure [5] shows how the various terms in the OPE con- 
tribute to the determination of the transition magnetic 
moment. The Af 2 term, which contains the contributions 
from the condensates x a an d strange quark mass m s , is 
the dominant term over the whole Borel region. The next 
leading contribution comes from Af term which involves 
condensates a and m s xa 2 , and is constant over the Borel 
region. The Af ~ 2 term is from condensates TOqO and %a&, 
and is slightly negative over the region. The Af ~ 4 term 
comes from condensates m s a 2 and ab, and is almost zero. 
These results confirm that the WO2 sum rule has good 
OPE convergence and why it is the best sum rule out of 
the three structures. It also shows the importance of the 
vacuum susceptibility x m the- transition. 

Figure[3]show various contributions in the phenomeno- 
logical representation in WO2 as a percentage of the to- 
tal. We need to make sure that the ground-state pole 
is dominant. We can see that it is about 70% at the 
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FIG. 2: Contributions from individual terms in the OPE for 
the WO2 sum rule. 



lower end of the Borel window. The transition contribu- 
tion is small in this sum rule. It is consistently smaller 
than the excited-state contribution and has a weak de- 
pendence on the Borel mass. The excited-state contribu- 
tion starts small, then grows with M 2 as expected from 
the continuum model. We cut off the Borel window when 
excited-state contribution goes too big, here it is around 
60%. 



70 - 




CT10- 



ts [ : 

1.6 1.8 2.0 2.2 2.4 2.6 2.8 

M 2 (GeV 2 ) 

FIG. 3: The three terms (pole, transition, and ESC) as a 
percentage of the spectral representation for the WO2 sum 
rule. 

To gain a deeper understanding of the underlying 
quark-gluon dynamics, it is useful to consider the indi- 
vidual quark sector contributions to the transition mag- 
netic moment. This can be done by dialing the corre- 
sponding charge factors e q in the sum rules which were 
explicitly kept for this purpose. Figure 0] shows a de- 
tailed matching of the two sides on Eq. (|24l) for the WO2 
sum rule, with individual contributions from u-, d- and 
s-quark. The four slopes are 0.80 for total fis°A an d 
0.53, 0.27 and for u-, d- and s-quark. We see that the 
u-quark contributes about 67% to the total transition 
magnetic moment, d-quark about 33%, and s-quark al- 
most zero. The slope for WO2 comes from 2 °£ , so 
it gives /is A — 1-60 fi^. To get the effective individual 
quark transition magnetic moment, we can use the per- 
centage of individual contribution from the total transi- 



tion magnetic moment. Then we have /i™a — 1-05 fiN, 
^ifiK ~ 0-53 hn and H^p^ ~ /ijv- A closer examination 
reveals that the smallness of the s-quark contribution is 
due to an almost exact cancellation of the terms involv- 
ing the s-quark, even though these terms are not small 
themselves. 
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FIG. 4: Individual quark contributions in the WO2 sum rule. 

We analyzed sum rules at structures WEi and WCq 
in a similar fashion. The full results determined at the 
structures are given in Table Q] We see that our calcu- 
lated transition magnetic moment is 1.60 ±0.07 /xjv from 
the most reliable sum rule W02- The other two sum 
rules also give close results around 1.60 /xjy. Our errors 
are derived from Monte-Carlo distributions which give 
more realistic estimation of the uncertainties. We find 
about 10% accuracy for the transition magnetic moment 
in our Monte-Carlo analysis, resulting from 10% uniform 
uncertainty in all the QCD input parameters. The uncer- 
tainties in the QCD parameters can be non-uniform. For 
example, we tried the uncertainty assignments (which are 
quite conservative) in Ref. fl3l | , and found about 30% un- 
certainties in our output. One choice for j3 here is —0.2 
because it minimizes the perturbative term in the mass 
sum rule [ljj ■ The excited state contribution at structure 
WO2 is 20-60%, which gives reasonable pole dominance 
over the whole Borel region. On the other hand, the sum 
rules at WEi and WOi structures have larger contribu- 
tions from the excited states over the Borel region. For 
this reason, we say the sum rules from WEi and WOi are 
less reliable. Another measure of reliability is the param- 
eter A, which is contamination from transitions between 
intermediate states. The A for the WOi sum rule is larger 
than those for the WEi and WO2 sum rules, making it 
less reliable. 



B. Relations between QCD sum rules 

It was pointed in Ref. [l8| that it is possible to derive 
QCD rum rules for A — > E° from those for S° and A. 
Here we would like to examine the issue using the gen- 
eralized interpolating fields for A and E° in Eq. ([3]) and 
Eq. Q. 
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TABLE I: Results for A — > £ from the three QCD sum rules. The ten columns correspond to, from left to right: structure, 
j3 value, Borel region in which the two sides of the QCD sum rule are matched, continuum threshold, percentage of excited 
state contribution on the phenomenological representation, transition strength, extracted transition magnetic moment in unit 
of nuclear magnetons, and effective u-, d- and s-quark transition magnetic moment. The errors are derived from 2000 samples 
in the Monte-Carlo analysis with 10% uncertainty on all QCD input parameters. 



Structure 


P 


Region 


w 


ESC 


A 




^£<>A 


,, d 


/4°A 






(GeV) 


(GeV) 


(%) 


(GeV- 2 ) 


(Mw) 


(Mn) 


(m«) 




WEI 


-0.2 


1.30 to 1.80 


1.80 


40-80 


0.11(1) 


1.57(8) 


1.04(8) 


0.53(3) 


0(0) 


WOl 


-0.2 


1.40 to 2.00 


1.70 


70-80 


-0.46(9) 


1.63(11) 


1.08(11) 


0.55(4) 


0(0) 


W02 


-0.2 


1.20 to 1.70 


1.70 


20-60 


0.17(4) 


1.60(7) 


1.05(10) 


0.53(3) 


0(0) 



The two-point correlation functions for magnetic 
moments of A and S° are constructed from time- 
ordered product of operators T{i] A (x) fj A (0)} and 
T{?7 S (x)ff 2 (0)}. The result is a master formula in 
terms of fully-interacting quark propagators in the QCD 
vacuum. Using the master formulas in Ref. Q, we have 
checked the following relations between the correlation 
functions 

2[nfl + n^j-n s ° = 3ir\ (25) 
2[n^ + n^j-n A = 3n s °. (26) 

Here fj^ s means the correlation function of A but with 
d-quark and s-quark interchanged. The relations imply 
that one can obtain the QCD sum rules for S° by sim- 
ple substitutions from the QCD sum rules for A and vice 
versa. We confirmed that relations not only at the corre- 
lation function level, but at the level of the derived sum 
rules in [8j. It serves as an independent check of the 
QCD sum rules in Ref. 

The same argument can be extended to the A — >■ S° 
transition moment. There exist the following relations 

jjs°(«« s ) _ -rj£ (<i«s) _ y3n s ° A ; (27) 

This suggests that one can derive the QCD sum rules 
for the A — > S° transition magnetic moment by simple 
substitutions, starting from those for either E° or A. In 
fact, we have calculated the QCD sum rules for A — > E° 
separately, and used the relations as independent checks 
of our calculations in this work and those in Ref. Q. 
This was done both at the level of the master formula in 
Eq. (fUtj) and at the level of the coefficients in the final 
QCD sum rules. The calculations has also served as a 
check of the Mathematica package (Tlj . 

C. Comparison of results 

Finally, we compare in Table |TT] our result with ex- 
periment and other calculations. The list is only a rep- 
resentative sample, and necessarily incomplete. Unlike 
the other members of the octet, there was only one 



TABLE II: Comparison of results for the A — > E° transition 
magnetic moments in nuclear magnetons. 



Method 


Me°a 


Experiment [19] 


1.61(8) 


QCD sum rule (this work) 


1.60(7) 


QCD sum rule [9] 


1.5 


Light-cone QCD sum rule [10J 


1.60 


Simple quark model [19, JO] 


1.57(1) 


Relativized quark model [21] 


1.04 


Quark model fits [22J 


1.81 


Chiral quark model [23] 


1.71, 1.66, 1.61 


Lattice QCD [24] 


1.36 


Chiral Perturbation Theory [25, 26] 


1.42(1), 1.61(1) 


Large-N ChPT [27] 


1.576 



measurement for the A — > E° transition magnetic mo- 
ment. It was done at Fermilab by utilizing the Pri- 
makoff effect [28j . Our result agrees with experiment. 
In the simple quark model, the transition moment can 
be related to the magnetic moment of E + and by 
A'soa = (/-*£+ — A i s-)v / 3/4. In terms of effective quark 
moments it is given by /is a = {p<u — (J-d)/V3. Using the 
values [i u = 1.852 /in and /i^ = —0.972 /ijv, it means 
that the u-quark contribution is about 67% of the total, 
the d-quark 33%, and the s-quark zero. Remarkably, our 
QCD sum rule results from the study of individual quark 
contributions are consistent with this composition. The 
lattice QCD result in Ref. [24| is quoted here by a scale 
factor of 1.18 as suggested in Ref. [29j |. 



V. CONCLUSION 

We have carried out a comprehensive study of the A 
— > S° transition magnetic moment using the method of 
QCD sum rules. We derived a complete set of QCD sum 
rules using generalized interpolating fields and analyzed 
them by a Monte-Carlo analysis. We determined the 
transition magnetic moment from the slope of straight 
lines. The linearity displayed from the OPE side matches 
almost perfectly with the phenomenological side in all 



9 



cases. Out of the three independent structures, we find 
that the sum rule from the W0 2 structure is the most re- 
liable based on OPE convergence and ground-state pole 
dominance. Our best prediction for the transition mag- 
netic moment is /i£°A = 1.60 ± 0.07 /j,n- More detailed 
results are listed in Table U The extracted transition 
magnetic moment is in good agreement experiment. Our 
Monte-Carlo analysis reveals that there is an uncertainty 
on the level of 10% in the transition magnetic moment 
if we assign 10% uncertainty in the QCD input param- 
eters. We find that for all three sum rules, the u-quark 
contributes about 2/3 of the total transition magnetic 
moment, d-quark about 1/3, while s-quark negligible due 
to an almost exact cancellation. This result agrees with 



the simple quark model. We verified certain relations 
that exist between the QCD sum rules and used them to 
check our calculations. Contrasted with results from a 
variety of other calculations, this study adds an unique 
perspective on the inner workings of the A — > S° transi- 
tion magnetic moment in terms of the non-perturbative 
nature of the QCD vacuum. 
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